Adaptive FEM with quasi-optimal overall cost for nonsymmetric linear elliptic PDEs
نویسندگان
چکیده
We consider a general nonsymmetric second-order linear elliptic PDE in the framework of Lax-Milgram lemma. formulate and analyze an adaptive finite element algorithm with arbitrary polynomial degree that steers mesh-refinement inexact iterative solution arising systems. More precisely, solver employs, as outer loop, so-called Zarantonello iteration to symmetrize system and, inner uniformly contractive algebraic solver, e.g., optimally preconditioned conjugate gradient method or optimal geometric multigrid algorithm. prove proposed iteratively symmetrized (AISFEM) leads full convergence for sufficiently small adaptivity parameters, rates respect overall computational cost, i.e., total time. Numerical experiments underline theory.
منابع مشابه
Adaptive FEM with Optimal Convergence Rates for a Certain Class of Nonsymmetric and Possibly Nonlinear Problems
We analyze adaptive mesh-refining algorithms for conforming finite element discretizations of certain non-linear second-order partial differential equations. We allow continuous polynomials of arbitrary, but fixed polynomial order. The adaptivity is driven by the residual error estimator. We prove convergence even with optimal algebraic convergence rates. In particular, our analysis covers gene...
متن کاملA heterogeneous stochastic FEM framework for elliptic PDEs
We introduce a new concept of sparsity for the stochastic elliptic operator −div (a(x,ω)∇(·)), which reflects the compactness of its inverse operator in the stochastic direction and allows for spatially heterogeneous stochastic structure. This new concept of sparsity motivates a heterogeneous stochastic finite element method (HSFEM) framework for linear elliptic equations, which discretizes the...
متن کاملSolutions to Problems for Quasi-Linear PDEs
(a) At what time ts and position xs does a shock first form? (b) Identify the important values of s and find the corresponding characteristics. Sketch the characteristics and indicate the region in the xt-plane in which the solution is well-defined (i.e. does not break down). (c) Construct tables as in class for the x and u values at the important values of s for times t = 0, 2/3, 4/3. Use thes...
متن کاملNew variants of the global Krylov type methods for linear systems with multiple right-hand sides arising in elliptic PDEs
In this paper, we present new variants of global bi-conjugate gradient (Gl-BiCG) and global bi-conjugate residual (Gl-BiCR) methods for solving nonsymmetric linear systems with multiple right-hand sides. These methods are based on global oblique projections of the initial residual onto a matrix Krylov subspace. It is shown that these new algorithms converge faster and more smoothly than the Gl-...
متن کاملEfficient Chebyshev Spectral Method for Solving Linear Elliptic PDEs Using Quasi-Inverse Technique
We present a systematic and efficient Chebyshev spectral method using quasiinverse technique to directly solve the second order equation with the homogeneous Robin boundary conditions and the fourth order equation with the first and second boundary conditions. The key to the efficiency of the method is to multiply quasiinverse matrix on both sides of discrete systems, which leads to band struct...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Ima Journal of Numerical Analysis
سال: 2023
ISSN: ['1464-3642', '0272-4979']
DOI: https://doi.org/10.1093/imanum/drad039